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Discovery of very- high- energy 7-ray emission from the vicinity of PSR J1831-952 
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We report on the latest discovery of an extended Very High Energy (VHE) 7-ray source near the 
67 ms pulsar PSR J1831-0952 during the H.E.S.S. Galactic Plane Survey (GPS). The dispersion 
measure distance of the pulsar (4.3 kpc) would imply that less than ~ 1% of its spin-down energy 
is required to provide the observed VHE luminosity of the source. No other plausible counterparts 
have yet been found through preliminary multi-wavelength searches. The most likely scenario is 
that the VHE emission originates from the - yet unseen at other wavelengths - wind nebula of PSR 
J1831-0952. If so this would constitute another case of a 7-ray discovered pulsar wind nebula. 



I. OBSERVATIONS AND ANALYSIS 

The H.E.S.S. (High Energy Stereoscopic System, 
an array of four imaging atmospheric Cherenkov tele- 
scopes located in the Khomas Highland in Namibia) 
has revealed more than 60 sources of very-high-energy 
(VHE) 7-rays through its Galactic Plane Survey 
(GPS) since 2004. Thanks to the use of advanced mul- 
tivariate analysis techniques and the accumulation of 
exposure, H.E.S.S. has achieved a sensitivity of better 
than 2% of Crab in the core region of the GPS (i.e. 
I = 282° to 60°) p. Pulsar Wind Nebulae (PWNe) 
constitute by far the dominating source population 
as compared to that of young shell-type Supernova 
Remnants (SNRs), or to that of older and/or inter- 
acting remnants. About one third of H.E.S.S. sources 
have still cither no known counterparts in other wave- 
lengths, or lack any clear emission scenario. 
The data on HESS J1831-098 consist of observations 
either dedicated to nearby sources such as SNR 21.5- 
0.9/HESS J1833-105, or from the extension of H.E.S.S 
GPS near I = 21°. These data were taken in 2004 
(May-Oct.), 2005 (June and July), 2007 (Apr. and 
July), 2008 (Sept.) and 2009 (Apr.-July), for a total 
observation time on HESS J1831-098 of ~ 52 hours. 
After application of the H.E.S.S. standard data qual- 
ity criteria [2] based on hardware and weather condi- 
tions, the data set for HESS J1831-098 amounts to a 
total live-time of 40 hours with average zenith angle 
and average offset (to the FoV center) of 22.8°, and 
1.30°, respectively. The mean offset is rather large 
because observations were not specifically targeted at 
this source. 



The standard Hillas H.E.S.S. event reconstruction 
scheme was applied to the data after calibration [5] ■ 
The rejection of cosmic-ray showers was done by ap- 
plication of a recently developed multivariate analysis 
The sky maps were produced with an image size 
cut of 80 photo-electrons (p.e.) and using the Ring 
Background method [5] where the background at each 
test position on the sky is derived from a ring sur- 
rounding it with a mean radius of 0.7°. To derive 
the spectrum, the same cut on image size was applied 
together with the a Reflected- Region procedure to es- 
timate the Background, followed by the application of 
a forward- folding method[3]. 



II. RESULTS 

The excess count map of the 0.4° x 0.4° region 
around HESS J1831-098 is shown on Figure 1. The 
map is smoothed with a Gaussian of a ~ 0.12°. An ex- 
tended 7-ray emission to the south-east of PSR J1831- 
0952 is observed with a peak pre-trials significance of 
7.9cr for the standard integration radius of 8 = 0.22°, 
used for generation of the GPS maps when search- 
ing for extended sources. The significance level of 
the source after a conservative correction for trials is 
5.8cr. The fit of the excess map with a two dimen- 
sional symmetrical Gaussian function, convolved with 
the H.E.S.S Point-Spread Function (PSF), results in 
a source centroid position of a ~ 18 h 31 m 25 s , 5 ~ 
-9° 54', with a width of a ~ 0.15° and a X 2 of 
593.4/525. The fit of an asymmetrical Gaussian func- 
tion does not improve significantly the % 2 nor the 
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FIG. 1: The excess map for HESS J1831-098, smoothed 
with a Gaussian of a ~ 0.12°. The white curves show the 
significance contours at 5, 6 and 7 a for an integration 
radius of 0.22°. The white cross shows the fitted position 
of the source. The position of PSR J1831-0952 is shown as 
a black triangle. Neighbouring SNRs are shown in green. 
Note that the fitted centroid does not coincide with the 
emission peak due to its departure from a Gaussian shape. 



FIG. 2: Differential energy spectrum of HESS J1831-098, 
fitted with a power-law function after extraction of events 
from a circular region of 0.3° radius, centered on the best- 
fitted position. The arrow shows the Fermi differential 
U.L. at 30 GeV. 



III. SEARCH FOR COUNTERPARTS 



residuals. A circular region of 0.3° radius, chosen 
as a compromise between optimal signal-to-noise ra- 
tio and independence of source morphology, was used 
to determine the energy spectrum. The differential 
spectrum, based on the 484 excess events in the circu- 
lar region, is shown in Fig. 2. The best fitted shape 
is a power-law (d</>/dE = </>o(E/lTeV)~ r ) with no sig- 
nificant indication of a cut-off up within the fitting 
energy range which extends up to 30 TeV. The fitted 
flux (uncorrected for events falling outside the integra- 
tion disk) is (j} = (1.1 ± 0.1) x lO-^TeV^cm^s- 1 
and the photon index T = 2.1 ± 0.1. The integrated 
flux for E > 1 TeV corresponds to about 4% of the flux 
of the Crab nebula in the same energy range [3J. 



A preliminary multi-wavelength search for energetic 
counterparts to HESS J1831-098 resulted in the sole 
pulsar PSR J1831-0952 [7 , which lies at a small an- 
gular offset of ~ 0.05° from the H.E.S.S. source's best 
fit position (see Fig. 1). Some X-ray data from Chan- 
dra and XMM- Newton are available, but the observa- 
tions have been done at large offsets with respect to 
PSR J1831-0952 and hence are not very useful. The 
search in the GeV energy was carried out using ~ 33 
months of public data from the Large Area Telescope 
(LAT) on board of Fermi (collected from 2008 Au- 
gust 4 (MJD 54682) to 2011 April 10 (MJD 55661)). 
A region of interest (ROI) of 6° was used to select 
events in the [10-100] GeV range. The ROI radius 
was chosen such as to be large enough to get a reli- 
able value for the normalization of the diffuse model 
and to be several times greater than the LAT PSF in 
the selected energy range. Events were analysed by 
applying the standard Fermi Science Tools to events 
of class 4 which is recommended for studies of faint 
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diffuse sources and which go beyond 20 GeV (in order 
to to minimize the non-photon background contami- 
nation) . The corresponding instrument response func- 
tion used for this event class is P6_V3_DATACLEAN 
[5]. Other standard cuts were also applied (e.g. zenith 
angle larger than 105° in order to reduce the contri- 
bution from terrestrial albedo 7-rays [3]). 

A hot-spot, at a pre-trials TS=25, was found within 
the boundaries of the H.E.S.S. source on a map gen- 
erated by gttsmap (Fig. 4). The hot-spot consists 
of one O.Oldeg 2 pixel at TS=25 and few neighbour- 
ing pixels at lower TS values. To account for the 
trials factor a conservative approach was adopted, 
namely the ratio of the predefined HESS source circle 
of radius 0.15° to one Fermi sky map pixel area was 
used. This yields a post-trial significance of ~ 4.5c 
for the hot-spot. Investigating a slightly higher en- 
ergy band of [30-100] GeV, the hot-spot significance 
drops to less than 3.4<r. Hence an upper limit (U.L.) 
on flux of U.L. ~ 5 x 10 _11 phcm _2 s _1 was derived by 
assuming a spectral index equal to that of the HESS 
source (r = 2.1). This U.L. is higher than the ex- 
trapolated flux of the H.E.S.S. source into the LAT 
range: ([30 - 100] GeV) - 3.5 x lO^phcm^s- 1 , 
and therefore does not exclude the existence of a GeV 
counterpart. 

IV. DISCUSSION AND SUMMARY 



FIG. 3: Fermi TS Map in the 10-100 GeV range (FoV = 
2°) with contours for TS > 30 corresponding to ~ 5.5a 
(pre-trials). HESS J1831-098 is shown by the yellow circle 
of radius 0.15°. 



PSR J1831-0952 is an energetic pulsar exhibiting a 
spin-down luminosity of 1.1 x 10 ergs , spin-down 
age of r c ~ 128 kyr, spin period of ~ 67 ms and a dis- 
tance estimated from dispersion measurements of 4.32 
kpc [10] . To power the H.E.S.S. source, a conversion 
efficiency from rotational energy to 1-20 TeV 7-rays 
of e ~ 1% would be required, i.e. a value compa- 
rable to those inferred for other VHE PWN candi- 
dates such as HESS J1420-609/HESSJ1418-607 in the 
wings of Kookaburra QI], HESS J1718-385 or HESS 
J1809-193 [12]. The angular offset of - 0.05°, trans- 
lating into a projected distance of ~ 4( jjr-j) pc, is also 
well within the range of offsets observed for crushed 
PWNe, as is the projected size of HESS J1831-098 
(~ 20( 4]^:) pc) when compared to the size of known 
VHE 7-ray PWNe. 

These values and the fact that the offset-type mor- 
phology is rather common to VHE emitting PWNe 
(e.g. HESS J1825-137, MSH 15-52, HESS 1718-385 
and HESS J1809-193) favour an interpretation of the 
VHE emission as originating from a PWN associated 
to PSR J1831-0952. In this scenario 7-rays are pro- 
duced through Inverse Compton (IC) scattering of 
ambient radiation fields (2.7 K Cosmic Microwave 
Background Radiation (CMBR), dust and star-light) 
by electrons injected by the pulsar into the nebula and 
re-accelerated at the wind terminal shock. The offset 
of the VHE-peak from the pulsar position can be due 



to the expansion of the SNR/PWN into an inhomoge- 
neous medium (e.g. [T3]) and/or the proper motion of 
the pulsar. In the latter case, and if the system age is 
indeed equal to the characteristic age of 128 kyr, the 
implied pulsar projected velocity of ~ 300 km s -1 re- 
mains reasonable and well within the bimodal velocity 
distribution derived in [T4"] . 

On the other hand, the extension of the maximum 
measured energy up to 30 TeV would imply a rather 
low magnetic field of 1 /iG in the IC scenario (see 
e.g., Eq. (6) in 15] ). However, the inferred value 
of the magnetic field intensity depends on the system 
age, which could suffer from large uncertainties (over- 
estimation) given that the pulsar characteristic age 
is calculated assuming that the initial spin period is 
negligible and remarking this is not always the case 
(e.g. PSR J1400-6326, where the true age of 1-2 kyr, 
is much smaller than the characteristic age of 12.7 kyr 
PI)- 

To summarize, although no PWN has been detected 
so far at other wavebands, the interpretation of HESS 
J1831-098 in terms of a wind nebula remains the most 
likely scenario, given the spatial coincidence with the 
energetic PSR J1831-0952, the reasonable efficiency 
and offset, and the observed abundance of such PWN- 
type VHE sources. If the association is confirmed, 
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HESS J1831-098 would constitute a gamma-ray dis- 
covered PWN. If so and if the true age of the system 
is close to the characteristic age of the pulsar, this 
source would be among the oldest known TeV PWNe. 
However, given the uncertainties on the source mor- 
phology due to limited statistics, other scenarios (e.g. 
a SNR shell emission) can not be excluded. X-ray ob- 
servations an additional data in the VHE 7-ray band 
are necessary to better understand the origin of the 
VHE emission. 
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